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Abstract — The gain flatness of EDFA (Erbium Doped Fiber 
Amplifier) plays an important role for WDM optical application 
and all optical self-routed wavelength addressable networks. 
EDFA have biggest disadvantage in having different gain for 
different wavelength. The main purpose of this work is to obtain 
the gain non-uniformity for each channel in order to equalize the 
amplitude gain in a WDM (Wavelength Division Multiplexing) 
system. The system is simulated to achieve gain flatness of EDFA 
through different fiber parameters. The design parameters for 
an EDFA simulation perspective are investigated. Complex 
effects of EDFA are researched and simulated with an optical 
environment design software tool named OptiSystem. 

Index Terms — ASE, BER, EDFA, WDM. 


the basic of a 980nm pumped EDFA. Consider the following 
three level system. 
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I. Introduction 

The new light wave generation, with vastly improved 
capacity and cost, is based on the recent development of 
EDFA’s. Undersea systems were the early beneficiaries, as 
EDFA repeaters replaced expensive and unreliable electronic 
regenerators. The use of Erbium doped fiber amplifiers is 
desirable in such systems to reduce the number of repeaters. 
An equally attractive feature of the EDFA is its wide gain 
bandwidth. Along with providing a gain at 1550nm, in the 
low-loss window of silicon fiber, it can provide a gain over a 
band that is more than 4000 GHz wide. With available WDM 
(Wavelength Division Multiplexing) technics and devices, 
commercial systems transport more than 16 channels on a 
single fiber; and the number is expected to reach 100. Thus, 
more than 250 optical channels with a channel spacing of 15 
GHz (0.1 nm) can be multiplexed into a fiber communication 
link and amplified by optical amplifiers placed periodically. 
But the EDFA gain varies with wavelength of signals. There 
comes the need of flattened EDFA gain and an optimized 
noise figure. 

In this study, we initially investigates the design parameters 
for an EDFA (Erbium Doped Fiber Amplifier) simulation 
perspective. A set of rate equations with boundary conditions 
are solved for the pump power inside the resonator cavity and 
the signal flow build up inside the cavity. Secondly, complex 
effects which occur during the gain and noise figure 
improvement operations of EDFA are researched and 
simulated with an optical environment design software tool 
named OptiSystem. 


II. Modelling 

For modelling an EDFA, we use the rate equations and 
propagation equations of a three level laser system which is 
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Figure 1. Three level system of laser 


Figure shows all the important transition between the levels. 
R- Absorption rate (pumping rate) from level 1 to level 3, 
corresponding to 980nm pumping. W i2 and W 2 i - Absorption 
rate and Stimulated emission rate between level 1 &2, A 21 - 
Spontaneous decay rate from level 2 to level 1, A 32 - 
Non-radiative decay rate from level 3 to level 2. By definition, 
level 1 is the ground level, level 2 is the metastable state 
having a life time x =1/ A 21 and level 3 is the pump level. Then 
the atomic laser rate equations of these three levels can be 
written as, 
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N 2 and jY 3 are the electron densities in levels 1,2 &3. By 
considering the fact that A 32 »R populations of the three 
level systems are, 
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Where, p = N 2 + jV 2 + A f 3 the erbium ion concentration. As 


shown in equation A r 3 is zero, this is due to the fast decay 
from level 3 to 2. Therefore 980nm pumped EDFA can be 
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analyzed by considering only levels 1&2. For a single mode 
fiber, R,W i2 and W 2 i, as a function of fiber coordinate z and 

radial axis r, can be written as, 
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tjj -denotes radial mode envelop distribution and subscripts s 
and p stands for pump and signal. o a , a e are the absorption and 
emission cross sections. Substituting eqns 7-9 in 5, we get the 
population density of level 2, 


Pump Power 



Figure 2. Pump power v/s fiber length 
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p(r) is the erbium density distribution. 


A. Propagation equation for signal 


When a signal of wavelength X with intensity I s (power per 
area) passes through an active medium of length dz, Then 
intensity change dl s , 

dh=(o s (X)N 2 - (JL)A \)I s dx (11) 

For a single mode propagation, and signal power of P s (^-), the 
light intensity distribution I s (A,r,0) 


tyXXr,9) 

I s (lr,d) - P S W Jiis(A ] riS y drdS (12) 

Where 0 is the azimuthal coordinate and S denote the integral 
should be over the entire transverse plane. From the above 
two equations, signal propagated is written as, 
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is the normalized mode power. 
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Above mentioned differential equations are solved by using 
fourth order Runge-Kuttae method to obtain signal and pump 
power flow in EDF. 


As fiber length increases, power gets reduced. 


Outut Power vs Pump Power 



Figure 3. Output power v/s pump power 

Graph shows pump and signal powers as a function of 
position along the fiber. Pump is injected at z=0 where the 
mirror reflectivity is high. 


Comparison: Pump and Sigal Powers 



III. Simulation 

The optical conversion efficiency is quite high (slope 
efficiency of 0.61) whereas the ratio of the output power to the 
absorbed pump power is 0.69. 


Figure 4. Backward signal inbuilt 

Pump and signal power as a function of position along the 
fiber, except that the pump is injected at z=L where the mirror 
reflectivity for the signal is much lower. 
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Comparison: Pump and Sigal Powers 




Figure 5. Forward signal inbuilt 

Although increasing the fiber length will absorb more power, 
at some point the signal output power will start to decrease 
due to the losses that become significant when the pump 
power is low. 


Figure 6. Use of gain flattening filter 

B. Changing Pump Power and EDF Length. 

Both pump power and EDF length is varied. Pump power is 
varied from 0 to 180mW. And EDF length is varied from 1 to 
6m. An optimized pump power and EDF length is found for 
the flattened gain EDFA. 


IV. Gain Flatening Methods Using Matlab 

In DWDM transmission systems and their related optical 
networks, one of the key technological issues is the 
achievement of broad and flat gain bandwidth for EDFA. 
Gain differences occur between optical channels having large 
wavelength spacing. In long amplifier chains, even small 
spectral gain variations can result in large difference in the 
received signal power, causing inacceptably large BER 
descrepancies between received signals. For some optical 
channels, complete power extinction can occur at the system 
output due to inefficient gain compensation along the 
amplifier chain. Additionally, the ASE generated in the region 
of highest gain in equalized EDF As causes homogeneous gain 
saturation, which affects WDM channels at longer 
wavelength. 



Dual Port WDM Analyzer | 

Lower frequency hoit = 185r TKz- 
Upper frequency limit = 200 THz 



Optical Power Meter Visuaiizer 
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WDM Transmitter Ideal Mux 

Frequency = 1558 nm 
Frequency spacing = -0.8 nm 
Power = -23.5 dBm 
Bit rate = Bit rate Bits s 
Modulation type = NRZ 




Optical Spectrum AnalyzerJ 


EDF 

Length = 4.754933456393365 m 


Frequency = 980 nm 

Power = 21.91829364728485 mW 
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Optical Spectrum Analyzer 


Figure 7. Changing pump power and EDF length 



Figure 5. Spectrum before flattening 
A. Use of Gain Flattening Filter. 

For a fixed length of EDFA, pump power is varied to achieve 
a predetermined gain. Output is optimized for different 
average insertion losses of Gain flattening filter. 


C. Adaptive method using matlab component 
Here we use a Matlab component as the gain flattening 
device. Code calculates the power obtained for each signal 
(wavelength). And selects the wavelength with minimum 
power.(or the wavelength got minimum gain).Then Matlab 
component attenuates other signal powers to the selected 
minimum power to equalize the gain. 
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m 


Optical Spectrum Analyzer_1 



Figure 6. Flattened spectrum 
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V. Conclusion 

In the design of EDFAs, many factors such as pumping 
wavelength, fiber length, fiber core radius, input signal 
power, fiber numerical aperture, and Erbium concentration 
have great influences. To find optimum values of these 
factors, we need a reliable and effective method. Although 
genetic algorithm can be quite helpful, sometimes it is hard to 
work with this method because of lengthy computational time 
and complexity. The population inversion can be controlled 
by optimizing fiber length and injected pump power to EDFA. 
With the proper selection of fiber length and injected pump 
power to EDFA, the population inversion can be controlled 
which in turn helps in maintaining the uniformity of gain in a 
WDM system. 
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